Cross Section Measurement of Charged Pion Photoproduction 
from Hydrogen and Deuterium 
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We have measured the differential cross section for the 771 — » -K~~p and jp — » ix + n reactions at 
0cm = 90° in the photon energy range from 1.1 to 5.5 GeV at Jefferson Lab (JLab). The data at 
E 1 > 3.3 GeV exhibit a global scaling behavior for both ir~ and tv + photoproduction, consistent 
with the constituent counting rule and the existing n + photoproduction data. Possible oscillations 
around the scaling value are suggested by these new data. The data show enhancement in the scaled 
cross section at a center-of-mass energy near 2.2 GeV. The cross section ratio of exclusive 7r~ to ir + 
photoproduction at high energy is consistent with the prediction based on one-hard-gluon-exchange 
diagrams. 



PACS numbers: 13.60.Le, 24.85. +p, 25.10.+S, 25.20.-x 

The study of the transition region from nucleon-meson 
degrees of freedom to quark-gluon degrees of freedom in 
quantum chromodynamics (QCD) for exclusive processes 
is one of the most interesting subjects in nuclear physics. 
Scaling in the differential cross section da/dt, and hadron 
helicity conservation have been pursued experimentally 
as signatures of this transition for years. While global 
scaling behavior has been observed in many exclusive 
processes, no experimental evidence supports hadron he- 
licity conservation. Furthermore, the exact nature gov- 
erning the onset of the scaling behavior is not clear. The 
relatively large cross section of pion photoproduction al- 
lows the search for additional possible signatures: QCD 
oscillations and the charged pion cross section ratio. In 
this experiment, three signatures (scaling, QCD oscilla- 
tions, charged pion ratio) for the transition are investi- 



gated. 

For an exclusive two-body reaction AB — > CD at high 
energy and large momentum transfer, the constituent 
counting rule (or the dimensional scaling law) predicts Q 

^-(AB^CD)~ s 2 - n f(9 cm ) , (1) 

where s and t are the Mandelstam variables. The quan- 
tity n is the total number of interacting elementary fields 
in the reaction, and f(9 cm ) is the angular dependence of 
the differential cross section. This rule was originally de- 
rived from dimensional analysis Q], and later confirmed 
within the framework of a perturbative QCD (pQCD) 
analysis up to a logarithmic factor of the strong cou- 
pling constant a s |2j. Many exclusive measurements at 
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fixed center-of-mass angles agree remarkably with this 

rule ma inn. 

The applicability of pQCD to exclusive processes re- 
mains controversial in the few GeV region. The pQCD 
calculation fails to correctly predict the magnitude of the 
proton magnetic form factor Furthermore, the onset 
of scaling can be sometimes as low as 1 GeV, as shown 
in the deuteron photodisintegration data 0, 0, which 
is much lower than the scale where pQCD is expected 
to be valid. Moreover, hadron helicity conservation, an- 
other consequence of pQCD (This statement is currently 
under debate 0), tends not to agree with polarization 
measurements at JLab in the photodisintegration pro- 
cess 0, E3, as well as in 1 HtfJ>)Tr , ed -> e d G3 
and ep — ► ep elastic scattering |lj| . 

In addition to the above puzzles, two striking anoma- 
lies have been observed in pp elastic scattering. The ra- 
tio of (da/dt)if /(da/ dt) n with spin normal to the scat- 
tering plane can reach 4 at 9 cm — 90° [j^j: the differ- 
ential cross section da/dt oscillates around the scaling 
value [la] • The interference between the short-distance 
and long-distance (Landshoff) subprocesses due to soft 
gluon radiation |l6j | can explain the above spin-spin cor- 
relation and oscillatory scaling behavior |l7j, |l8( . This 
interference is analogous to the QED effect of coulomb- 
nuclear interference observed in charged particle scatter- 
ing at low energy. Alternatively, the above anomalies in 
pp scattering can be interpreted in terms of resonances 
associated with charm production threshold, interfering 
with a pQCD background [lflj . 

It was previously thought that the oscillatory scaling 
behavior was unique to pp scattering or hadron induced 
exclusive processes. However, it has been suggested 
that similar oscillations should occur in deuteron pho- 
todisintegration [20| and pion photoproduction at large 
center-of-mass angles |21|. The recent photodisintegra- 
tion d(-f,p)n data 0, 0| showed that the oscillations if 
present are very weak, and the rapid decrease in the cross 
section with photon energy (da/dt oc s~ n ) makes it im- 
practical to investigate such oscillatory behavior. Thus, 
it is essential to search for oscillations in pion photopro- 
duction, which has a much larger cross section at high 
energy due to a slower decrease in the cross section with 
energy (da/dt oc s~ 7 ), compared to the deuteron pho- 
todisintegration. In this Letter, we present cross section 
results for charged pion photoproduction from hydrogen 
and deuterium at 9 cm = 90°. This angle was chosen to 
achieve the highest transverse momentum, which might 
be the kinematic quantity indicative for the onset of seal- 
ing 0,113. 

The reaction in which pQCD is expected to manifest 
itself at relatively low momentum transfer is the elas- 
tic scattering from the charged pion, since the pion has 
the simplest valence quark structure. The study of the 
charged pion elastic form factor has thus been pursued 
vigorously in recent years in both experiment and the- 



ory [U 0| . The elastic form factor of the charged 
pion is still dominated by non-perturbative effects at a 
few (GeV/c) 2 . But these non-perturbative effects cancel 
out to first order in the differential cross section ratio 
( 7TW7r p ) Therefore one may expect the ratio to 

give the first indication of a simple pQCD prediction. In 
this letter, we provide data on the charged pion ratio up 
to a momentum transfer squared value of 5.0 (GeV/c) 2 , 
the highest value ever achieved for this quantity. 

Experiment E94-104 was carried out in Hall A [2^] 
at the Thomas Jefferson National Accelerator Facility 
(JLab). The continuous electron beam, at currents 
around 30 fiA and energies from 1.1 to 5.6 GeV, impinged 
on a 6% copper radiator and generated an untagged 
bremsstrahlung photon beam. The production data were 
taken with the 15 cm cryogenic liquid hydrogen (LH2) 
target for singles 75(7, 7r + )n measurement, or with the 
liquid deuterium (LD2) target for coincidence ^(7, ir~p)p 
measurement. The background was measured by record- 
ing data without the production target and with or with- 
out the radiator. The two High Resolution Spectrome- 
ters (HRS) in Hall A, with a momentum resolution of 
better than 2 x 10~ 4 and a horizontal angular resolution 
of better than 2 mrad, were used to detect the outgoing 
pions and recoil protons. Based on two-body kinematics, 
the incident photon energy was reconstructed from final 
states, i.e. the momentum and angle of the tt + in the 
singles measurement, momenta and angles of the ir~ and 
p in the coincidence measurement. Both spectrometers 
consisted of magnets to focus and bend the charged par- 
ticles (45°), vertical drift chambers (VDCs) to record the 
tracks, and scintillator planes (S1/S2) to generate trig- 
gers. Two new aerogel Cerenkov detectors (A1/A2) in 
the left spectrometer provided particle identification for 
positive particles, mainly pions and protons, since the 
time-of-flight technique fails at high momentum. The 
CO2 gas Cerenkov detector and preshower/shower de- 
tector in the right spectrometer provided particle iden- 
tification for negative particles, mainly pions and elec- 
trons. The average number of detected photoelectrons 
generated by a (3 — 1 particle was approximately 8 for 
the aerogel Cerenkov detector Al (n — 1.015), approxi- 
mately 30 for A2 (n = 1.055) and approximately 10 for 
the gas Cerenkov detector (n = 1.00041). 

A 100 MeV bin of the reconstructed photon energy 
spectrum, centered 75 MeV below the beam energy, was 
chosen for the data analysis, where the multi-pion contri- 
bution is negligible. The data after background subtrac- 
tion, with cuts on trigger type, coincidence timing, PID 
(particle identification), acceptance and photon energy, 
were compared to the Monte Carlo simulations to extract 
the raw cross section. The simulation was done with a 
JLab Hall A Monte Carlo program, MCEEP ,26|, with 
some modifications. The bremsstrahlung photon flux was 
calculated with an estimated 3% uncertainty, by using 
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the thick-radiator codes written by Meekins |27j, based 
on the formulae of Matthews et al. [2^] . The momentum 
distribution of the neutrons in the deuterium target as 
well as the binding energy was considered in the simu- 
lation. Different calculations of the momentum distribu- 
tion were tried and little model dependence was found for 
the cross section (< 1%). The angular distribution input 
for the cross sections was fitted from the ir + photopro- 
duction data at 4, 5 and 7.5 GeV 0, which has not been 
reproduced by the pQCD calculation j^. It was used 
for all the kinematics including the 7r~ photoproduction 
since the extraction of cross section at 6 cm — 90° was 
insensitive to the angular distribution (< 1%). The dis- 
tributions of acceptance, reconstructed momentum and 
photon energy from data were in good agreement with 
those from simulations. 

Several correction factors were applied to deduce the 
final cross section, as shown in Table |l] The largest cor- 
rection, on the order of 20%, came from the the nuclear 
transparency in the deuteron due to final state interac- 
tions. The nuclear transparency was obtained for the 
d(j,ir~p)p process based on a Glauber calculation [30| . 
which has been tested by the measured transparency 
from the quasi-elastic d(e,e'p) process The correc- 
tion due to material absorption was applied to compen- 
sate for the scattering losses in the target and in the 
spectrometers. The correction for a single pion or proton 
was approximately 6%, with major losses in the target, 
scintillators and aerogel detectors. The pion decay losses 
were calculated from the flight distance. As some muons 
from pion decay may still fall into the acceptance and be 
mis-identified as pions, an additional correction has to be 
applied, which was 4~7% based on Monte Carlo simula- 
tions [23| . The computer dead time correction was con- 
sidered run- by-run and was mostly a few percent. The 
detector efficiencies also led to some corrections, mostly 
less than 1%. 

The total errors were dominated by systematic uncer- 
tainties, which were estimated to be 8% in cross section. 
The point-to-point systematic uncertainty for the three 
kinematics at 3.3, 4.2 and 5.5 GeV is 4%. The sta- 
tistical errors were approximately 2%. The major sys- 
tematic uncertainties arose from the calculation of the 
bremsstrahlung photon yield, the simulation of the ac- 
ceptance and the estimation of the nuclear transparency, 
material absorption and pion decay factor, approximately 
3% for each item. Also, there were 2% uncertainties from 
PID and the energy loss calculations. Other systematic 
uncertainties were less than 1%. 

The upper plot in Fig.^shows the results of the scaled 
differential cross section (s 7 ^ 1 ) for the jp — > n + n pro- 
cess at 9 cm = 90°. The new results with fitted value 
n = 9.0 ± 0.2 (see Eq. 1) agree with those of Anderson 
et al. 4] and exhibit the scaling behavior predicted by 
the constituent counting rule with 9 elementary fields. 
The lowest energy datum in the inset box of Fig.n cor_ 
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FIG. 1 

center-of-mass energy for the jp — > iv + n (upper plot) and 
7n — > iv~p (lower plot) at 6 cm = 90°. The data from JLab 
E94-104 are shown as solid circles. The error bars for the 
new data and Anderson et al.'s data Q, include statistical 
and systematic uncertainties, except that those in the insets 
only include point-to-point uncertainties to highlight the pos- 
sible oscillatory scaling behavior. Other data sets [33l . |34| 
are shown with only statistical errors. The open squares in 
the lower plot were averaged from data at 8 cm — 85° and 
95° |3 51) - The solid line was obtained from the recent partial- 
wave analysis of single-pion photoproduction data |3al up to 



E 7 =2 GeV, while dash line from the MAID analysis |37rj up 
to E^=1.25 GeV. 



responds to a center-of-mass energy of approximately 2.7 
GeV and photon energy of 3.3 GeV. The corresponding 
transverse momentum is approximately 1.2 GeV/c. 

The lower plot in Fig^shows the results of the scaled 
differential cross section (s 7 for the jn — * n~p pro- 
cess at 9 cm = 90°. The new results greatly extend the 
existing measurements and exhibit, for the first time, a 
global scaling behavior at high energy for this reaction 
with fitted value n = 8.6 ± 0.2. The scaling behavior 
in 7r~ production is similar to that in 7r + production. 
Furthermore, data in these two channels show possible 
oscillations around the scaling behavior in similar ways 
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as suggested by the insets in Fig. ^ Note that this possi- 
ble oscillatory behavior occurs above the known baryon 
resonance region. Unfortunately, the coarse and few pho- 
ton energy settings of this experiment do not allow us to 
claim the observation of oscillations. Measurements with 
much finer binning, planned at JLab (3^, are essential 
for the confirmation of such oscillatory scaling behavior. 

The future 12 GeV energy upgrade at JLab will en- 
able us to extend further the measurements and study 
the reactions both below and above the charm produc- 
tion threshold. Searching for the resonances around the 
charm production threshold, proposed by Brodsky and de 
Teramond 0] to explain the oscillatory scaling behavior 
in pp scattering, will help to understand the origin of the 
oscillatory scaling behavior. 



the exact nature of the observed enhancement. 

The cross section ratio of n~ to ir + photoproduction 
can be calculated 01 based on one-hard-gluon-exchange 
diagrams as 

dajjn -> n~p) _ ( ue d + se u \ 2 
dcr(jp^ir+n) \ue u + se d J 

where u and s are the Mandelstam variables, and e q de- 
notes the charge of the quark q. The non-perturbative 
components are represented by the form factors which 
divide out when the ratio is taken. The calculation is 
expected to be valid only at high energy. As shown in 
Fig. [3J the calculation agree with the two data points at 
the highest energies. 



E-y 

(GeV) 


f ( 7 p^7r+n) 
(nb/GeV 2 ) 


^( 7 n->7r p) 
(nb/GeV 2 ) 


1.106 


(1.16±0.01±0.09)xl0 4 


(5.72±0.03±0.46)xl0 3 


1.659 


(1.36±0.01±0.11)xl0 3 


(2.39±0.01±0.19)xl0 3 


1.815 


(1.06±0.01±0.08)xl0 3 


(1.58±0.01±0.13)xl0 3 


2.481 


(1.87±0.02±0.15)xl0 2 


(2.43±0.03±0.19)xl0 2 


3.321 


8.07±0.09±0.65 


(l.ieio.oiio.ogjxio 1 


4.158 


2.34±0.04±0.19 


4.05±0.08±0.32 


5.536 


0.33±0.02±0.03 


0.56±0.01±0.04 



TABLE I: The differential cross section ^ at 8 cm = 90° for 
7P — > TT + n and -yn — > ir~p reactions followed by the statistical 
and systematic errors. 



Another interesting feature of the data is an apparent 
enhancement in the scaled differential cross section be- 
low the scaling region, at a center-of-mass energy rang- 
ing approximately from 1.8 GeV to 2.5 GeV, in both 
channels of the charged pion photoproduction, as shown 
in Fig. ^ This effect was also observed in neutral pion 
photoproduction US ESI- Without any conclusive state- 
ments for the present, some speculations can be made. 
The observed enhancement around 2.2 GeV might relate 
to some unknown baryon resonances, as some of the well 
known baryon resonances (A, N*'s around 1.5 GeV and 
1.7 GeV) are clearly seen in the scaled cross section below 
2.2 GeV. Several baryon resonances are predicted to be 
in this energy region by the constituent quark model [13 j 
but have not been seen experimentally, i.e. the so called 
'missing resonances'. The observed enhancement might 
be associated with the strangeness production threshold. 
When the available energy is near a new flavor threshold, 
all of the quarks have small relative velocities, allowing 
resonant behavior in attractive channels Whatever 
the origin is for the observed enhancement around 2.2 
GeV, it is remarkable that the scaled differential cross 
section drops by a factor of 4 over a rather small center- 
of-mass energy range (~ 0.3 GeV). More experimental 
work, such as polarization measurements |4(J, together 
with an extended SAID model might help to understand 




FIG. 2: The cross section ratio of n to n photoproduction 
versus center-of-mass energy 

In summary, we have measured the differential cross 
section da/dt for the photoproduction processes of jn 
tt~P with a deuterium target and jp — > 7r + n with a hy- 
drogen target at 8 cm = 90° with photon energies from 
1.1 to 5.5 GeV. The data with E 1 > 3.3 GeV exhibit a 
global scaling behavior in both processes, consistent with 
the constituent counting rule. The data with E 7 > 3.3 
GeV also suggest a possible oscillatory scaling behav- 
ior, the confirmation of which awaits future measurement 
with finer binning in energy. Furthermore, the scaled 
cross section data show an enhancement at a center-of- 
mass energy near 2.2 GeV and the exact nature of such 
a structure requires further investigation. The data also 
provide ir~ to ir + cross section ratios, consistent with the 
one-hard-gluon-exchange prediction at high energies. 
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